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en ' Abstract 

I ' We derive explicit transformation formulae relating the renormalized 

r^ , quark mass and field Eis defined in the MS-scheme with the correspond- 

0\ ' ing quantities defined in any other scheme. By analytically computing 

the three-loop quark propagator in the high-energy limit (that is keeping 

only massless terms and terms of first order in the quark mass) we find 

the NNNLO conversion factors transforming the MS quark mass and the 

^ l' renormalized quark field to those defined in a "Regularization Invariant" 

D . (RI) scheme which is more suitable for lattice QCD calculations. The 

i-C ' NNNLO contribution in the mass conversion factor turns out to be large 

and comparable to the previous NNLO contribution at a scale of 2 GeV 
— the typical normalization scale employed in lattice simulations. Thus, 
in order to get a precise prediction for the MS masses of the light quarks 
from lattice calculations the latter should use a somewhat higher scale of 
around, say, 3 GeV where the (apparent) convergence of the perturbative 
series for the mass conversion factor is better. 

We also compute two more terms in the high-energy expansion of the 
MS renormalized quark propagator. The result is then used to discuss the 
uncertainty caused by the use of the high energy limit in determining the 
MS mass of the charmed quark. As a by-product of our calculations we 
determine the four-loop anomalous dimensions of the quark mass and field 
in the Regularization Invariant scheme. Finally, we discuss some physi- 
cal reasons lying behind the striking absence of ^(4) in these computed 
anomalous dimensions. 
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1 Introduction 

Quark masses are fundamental parameters of the QCD Lagrangian. Neverthe- 
less, their relation to measurable physical quantities is not direct: the masses 
depend on the renormalization scheme and, within a given one, on the renor- 
malization scale /i. 

In the realm of pQCD the definition which is most often used is based on 
the MS-scheme [|l|, D which leads to the so-called short-distance MS mass. Such 
a definition is of great convenience for dealing with mass-dependent inclusive 
physical observables dominated by short distances (for a review see H). Un- 
fortunately it is usually difficult to get precise information about the quark 
masses from predictions from these considerations, as their mass dependence is 
relatively weak. 

To determine the absolute values of quark masses, one necessarily has to 
rely on the methods which incorporate the features of nonperturbative QCD. 
So far, the only two methods which are based on QCD from first principles are 
QCD sum rules and Lattice QCD (for recent discussions see e.g. M H, h, u[ N, 0, 
00, O, 12, Of. Rather accurate determinations of the ratios of various quarks 



masses can be obtained within Chiral Perturbation Theory ll£ 

Lattice QCD provides a direct way to determine quark masses from first 
principles. Unlike QCD sum rules it does not require model assumptions. It 
is possible to carry out a systematic improvement of Lattice QCD so that all 
the discretization errors proportional to the lattice spacing are eliminated (a 
comprehensible review is given in [15|). The resulting quark mass is the (short 
distance) bare lattice quark mass. The matching of the lattice quark masses 
to those defined in a continuum perturbative scheme requires the calculation of 
the corresponding multiplicative renormalization constants. In the RI scheme 
[EgI the renormalization conditions are applied to amputated Green functions 
in Landau gauge, setting them equal to their tree-level values. This allows the 
non-perturbative calculation of the renormalization constants. An alternative 
to the RI approach is the Schroedinger functional scheme (SF) which was used 
in0,p. 

An impressive number of various lattice determinations of quark masses 
recently has been performed (see Refs [||, |^, |l[ || ||. ||, || ^. |^ ||, |§, 
013). 

Once the RI quark masses are determined from lattice calculations they can 
be related to the MS mass by a corresponding conversion factor. By neces- 
sity this factor can be defined and, hence, computed only perturbatively . The 
conversion factor is presently known at next-to-next-to- leading order (NNLO) 
from IQ. The NNLO contribution happens to be numerically significant. This 
makes mandatory to know the NNNLO 0{al) term in the conversion factor. 

In the present article we describe the calculation of this term. It turns out 
that the size of the newly computed term is comparable to the previous one at 
a renormalization scale of 2 GeV — the typical scale currently used in lattice 
calculations of the light quark masses. This means that perturbation theory can 
not be used for a precise conversion of the presently available RI quark masses to 



the MS ones. A simple analysis shows that the convergence gets much better if 
the scale is increased to, say, 3 GeV. Thus, once the lattice calculations produce 
the RI quark masses at this scale our formulas will allow an accurate conversion 
to the MS masses at the same scale. 

The article is organized as follows. In section we discuss the scheme de- 
pendence of the quark field and mass and present a general procedure to find 
corresponding conversion factors from one scheme to another. The general tech- 
nique is illustrated by constructing the conversion factors between the MS and 
the RI scheme. In section |3| we present the first few terms of the small mass 
expansion of the three-loop quark-propagator in the MS scheme. In section Q 
we first present the conversion functions for the quark mass and field, and then 
investigate the validity of the massless approximation for these functions. Then 
we use these results to calculate the anomalous dimensions for the quark mass 
and field in the RI scheme, the so-called RG invariant mass m. The final section 
is devoted to conclusions. 

In Appendices A and B we display our results for the small mass expansion 
of the fermion propagator and the various conversion factors with their full 
dependence on the group theoretical factors Cp, Ca and T. In Appendix C the 
four loop anomalous dimensions of the quark mass and field are listed for the 
case of a SU(A^) gauge group. 

Our main results are also available as ASCII input for the programs FORM 
and Mathematica at the following internet address: 



tittp : //www-ttp .physik.uni-karlsruhe .de/Progdata/ttp99/ttp99-43/ 



2 Scheme dependence of quark mass and field 

2.1 Generalities 

We start by considering the bare quark propagator (for simplicity we stick to 
the Landau gauge in this section and, thus, do not explicitly display the gauge 
dependence) 

S^{q,almo) = i Uxe^'i^i:[Mx)Mm = ^^^, (1) 

J mo-^-i^o 

with the quark mass operator Eq being conveniently decomposed into Lorentz 
invariant structures according to 

j:o = ^^v+mo^S^ (2) 

where mo and V'o are the bare quark mass and field, respectively. Additionally 
we are using the common shortcuts for the coupling constant throughout this 
article: 

_ «° ^ ffo 

go being the bare QCD gauge coupling. To be precise we assume that (nh is 
dimensionally regulated by going to non-integer values of the space-time di- 
mension Z) = 4 — 2e pS, B4[. The MS renormalizcd counterpart of the Green 



function (m reads 

S{q,as,m,fi) = i / da:e''-(T[V'(x)V;(0)]) = 



m, — g — E 
= Z^^So{q,a°,ma)\ o ,^ , (3) 

where the renormaUzed quark field is 

and the 't Hooft mass parameter [i, is the scale at which the renormalized quark 
mass is defined. The renormalization constants Z2, Z^ and Z^ are series of the 
generic form 

z.^i+^z<^)^ ,z«.5:z(-)(^y. (4) 

i>0 j>i 

A general theorem (first rigorously proven for minimal subtraction in ||35| , |36|) 
states that there is a unique choice of renormalization constants of the form (^ 
which makes the propagatoi|j finite in the limit of _D — > 4. 

The independence of the bare coupling constant, mass and quark field on /i 
leads in the standard way to the following renormalization group equations for 
their renormalized counterparts 



^ i>0 i>0 



Pias) ^ -J2P.{^y^" =Y.(^)Z^»'^^ (5) 



M^^"^(A')lgo,mo = m(^)7„(Q;,) = -771^7^ (^_ij 



i+l 



i>0 

(1) 



^ "^(a^)«.^ (6) 



and 



= i^i^^)as^. (7) 

Now let us consider the quark propagator renormalized according to a dif- 
ferent subtraction procedure. Marking with a prime parameters of the second 
scheme one can write 

S'{q,a',,m',^^) = i f dxe''^^T[i^'{x)i^'{Q)]) ^ — ^ 



m' — ^ — S' 



iZ!,)-'So{q,almo)\ „ , , (8) 



^ In fact all Green functions if proper renormalization constants for gluon and ghost fields 
are introduced. 



where without essential loss of generality we have set //' — fi. The finiteness of 
the renormalized fields and parameters in both schemes implies that, within the 
framework of perturbation theory, the relation between them can uniquely be 
described as follows 

fn — Cm ■ m' (9) 

V' = \/^-V'', (10) 

with the "conversion functions" being themselves finite series in aj,, i.e. 

c..i+5:c«(^)^ (11) 

for ? = TO or 2. 

Note that in general the coefficients C? may depend on the ratio m' / jjl. 
If such a dependence is absent then the corresponding subtraction scheme is 
referred to as a "mass independent" one. In what follows we mainly limit 
ourselves to considering this latter case. In addition, being only interested in 
the conversion functions C2 and Cm, we will assume that the function Ca has 
already been determined and, thus, will deal with the following representation 
of C2 and Cm in terms of the MS coupling constant as'. 

C,^1 + 5]C«(^)\ (12) 

The running of m' and ■0' is governed by the corresponding anomalous di- 
mensions 7m(as) and 72(as)- A direct use of Eqs. (||,|o|) gives 

d 

im=l^n- Pj— In Cm , (13) 

oaa 

i2=l2-Pj-\^C2. (14) 

das 
At last, from Eq. dl^) it is easy to see that 

^('^) = ^-^'(-^) = TO-(l^E-/!^(l + EV) ^''^ 

or, equivalently, 

C2 • (1 + Sy) = 1 + SV (16) 

C2 ■ Cm • (1 - Ss) = 1 - S's. (17) 

The renormalization conditions for the non-MS scheme should then be used 
to provide the necessary information about the right hand side to calculate the 
conversion factors Cm and C2 once the MS renormalized Ey and Eg are known. 



2.2 Regular izat ion Invariant scheme versus MS 

The MS subtraction scheme is intimately connected to dimensional regulariza- 
tion and, thus, can not be directly used with other regularizations, including the 
lattice one. In addition, the physical meaning of its normalization parameter 
H is not transparent and leads to the well-known ambiguities when considering 
the decoupling of heavy particles. 

It is well-known that the above shortcomings are absent for a wide class of 
so-called momentum subtraction (MOM) scheme^ The MOM schemes require 
the values of properly chosen Green functions with predefined /j, dependent 
configurations of external momenta to be fixed (usually to their tree values) 
independently on the considered order. Practical calculations can then be per- 
formed with any regulator (or even without it in the regulator-free approach 
of [^ |38|). A shortcoming of MOM schemes is that they are in general not 
mass-independent which leads to a complicated running of coupling constant (s) 
and mass(es). 

A general analysis of the problem of constructing mass-independent subtrac- 
tion schemes was performed long ago in [p^. Following essentially Weinberg's 
ideas, a specific example of a mass independent MOM scheme for QCD has 
recently been considered in [|l6| under the name of RI ( "Regularization Invari- 
ant") scheme (see also pH]). The corresponding renormalization conditions for 
i:f and S«;i read 



lim 4^Tr 



m^O 



48' 



lim TTjTr 



7p 



a(^(i + s^^)) 



1 ^^S\2=_^2 = 1, 



= 1 



(18) 



where the trace is to be taken over Dirac, Lorentz and colour indices. Note 
that the zero mass limit in (|lj) means that both Sy' and Sg' are effectively 
massless functions only depending on the QCD coupling constant, the normal- 
ization point /i and q^ . This also implies that it is sufficient to compute the MS 
functions Ey and Eg in massless QCD when computing the conversion factors 
from relations (16h and (p^). 

Application of the renormalization conditions (Eq) to the conversion formulae 
(Iq) and (O) leads to equations that can simply be solved for C^ and C^^ . All 
the dependence of Ey on q^ is of the form of £ = log{--^), which simplifies the 
trace and derivative w.r.t. g^ and leads to 



Cf = 



1 dJ^vji) 



(19) 



-^2^ 771 — 






RI 



1 d^vji) 

2 • d(! 



1-Es 



(20) 



-fi^ m— 



■^ In a sense the oldest subtraction scheme — the on-shell one for QED — can also be 
considered as an example of a MOM scheme. 



Another quark field renormalization that is useful in numerical lattice simula- 
tions is defined by the condition 



hm —Tr 

m^O 12 



1+sr 



1, 



which results in the even simpler conversion factors 

-/^^, m— 



cf = [i + Sv];.L 



G 



RI' 



1 



^V 



1 



(21) 



(22) 



(23) 



-/A^, m— 



Using these equations all conversion factors can easily be obtained, once the 
MS renormalized expressions Ss and Ey are known. 

It should be noted that in practical lattice calculations the massless limit 
on the left hand side of (|l8) and ( pl| ) is implemented by choosing fi ^ m. On 
the other hand, fi should be much less than the inverse lattice spacing 1/a. 
Typically /i is taken around 2 GeV. This means that lattice determinations do 
not lead directly to the RI quark mass but rather to the mass in a different, mass- 
dependent, scheme. The difference between both schemes can be numerically 
non-negligible in the case of the charmed quark. 

Thus, for both the RI and the RI' scheme it is suggestive to introduce their 
mass-dependent counterparts MOM and MOM' as defined by the same Eqs. ( [l8[) 
and ( pi[) but without taking the m ^ Q limit. The corresponding conversion 
factors to the MS scheme read 



a 



MOM 



a 



MOM 



1 25Sy 



l + Sv + i-q2^ 



and 



MOM' 



C2 

/^MOM' 



I-S5 

[1 + I]y]g2^_^2 
1 + Sy' 



1-S.' 



(24) 

(25) 

(26) 
(27) 



respectively. 



3 Three-loop quark propagator in MS-scheme 

To find the conversion factors for MOM and MOM' one needs to compute the 
functions Sy and Eg including their full mass dependence. A full analytical 
result at two-loop level has been obtained only recently in U^ . An extension of 



this calculation up to three-loops is out of reach of present calculational tech- 
nologies. Fortunately for the RI and RI' schemes one effectively only needs to 
compute massless three-loop diagrams - a problem which in principle was solved 
long ago in p2| . A promising approach to recover the full mass dependence of 
the quark propagator seems to be to employ an expansion in {rn? /q^) p3[ . In- 
deed, as has been demonstrated in H, E^, El, EtI small mass expansions can 
be a very effective tool for accurate predictions of mass dependences provided 



one is not too close to the threshold {q 



in our case). The exact two- loop 



result for the quark propagator can provide some insight into the accuracy of 
such an expansion. 

We have analytically computed three terms in the small mass expansion of 
the quark propagator to order a^. The calculation has been done with intensive 
use of computer algebra programs. In particular, we have used QGRAF |Q for 
the generation of diagrams and LMP [^ for the diagrammatic small mass ex- 
pansion. The small mass expansion results in products of massless propagators 
and massive tadpoles. These have been evaluated with the help of the FORM 
packages MATAD |^ and MINCER ^^ (A detailed description of the status of 
these algebraic programs can be found in [ p2| ). It is convenient to write the 
functions Ss/v ii^ the following way: 



^s/v = 1 






(i) 

S/V 



Um) 



(28) 



Our results for the separate contributions in the Landau gauge and in the MS 
scheme read, where we here only keep the rif dependence, n/ being the number 
of light quark flavours, with one quark flavour of mass m and Uf ~ \ massless 



quark flavours (Zq„ 
Landau gauge Sy : 



\n{-q^/m^), Iq 



0): 



ln(— (J //i ), z = m /q , and in 
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48'^'^ 
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216 
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- ^ ^ff'qm 



18 



Uf 
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'ou 1 



—li 
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g ^At^ 



109, 
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209 
144 



gauge parameter and the group 
Note that the result for the three 



•^the expressions including the full dependence on th 
theoretical factors Ca , Cp and T are given in Appendix [ 
loop massless quark propagator can also be found in the source code of the updated version 
of the FORM program MINCER jHsl. 
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4 RI scheme versus MS and MOM schemes 

4.1 Three Loop Conversion functions 

A direct use of Eqs. (|l^,^,^^ and (|23| ) leads to the following analytical expres- 
sions for the conversion factors between the MS and RI schemes. The results 
are shown for QCD (SU(3)) and Landau gauge as functions of n/ ( the Q are 
the values C(i) of Riemann's Zeta function ): 
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(37) 



At a scale oi fi ^ 2 GeV and nj — 4, the numerical contributions of the 
leading order to NNNLO terms are as follows (for simplicity we inserted Os/tt = 
0.1): 






RI 



RI 



1.0 + 0.0 - 0.00476 - 0.00508 , 
1.0 - 0.1333 - 0.0754 - 0.0495 



(38) 
(39) 



and 



Cf'' = 1.0 



C^'' = 1.0 



0.0-0.0101-0.0095, (40) 

0.1333-0.0701-0.0458. (41) 

One observes that the sizes of the NNLO and NNNLO contributions to C^ at 
this scale amount to about 7.5% and 5% respectively. This shows that pertur- 
bation theory can not be used for a precise conversion of the RI quark masses 
to the MS ones at the renormalization scale fj, = 2 GeV. The convergence can 
be improved if one increases fj, to, say, 3 GeV. Indeed, with this choice of ^ 
the standard three-loop evolution gives as (3 GeV) — 0.262 and Eqs. ( p8| , p9| ) 
transform to 



and 



Cf = 1.0 



0.0-0.00333-0.00296 



a" = 1.0 - 0.111 - 0.0526 - 0.0289. 



(42) 
(43) 
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The accuracy of the massless approximation can be tested by computing the 



ratio Cf^/a 



MOM 



{'■ 



, 2) as a series in z = miK—y?). Using the results of 



the previous section one obtains 



(^MOM 



[0.28981 z ~ 0.89236 z^ + 1.5284 z^ - 3.3649 z"* + 7.6945 z^ 
+ 0.25 zh - 0.39583 z^l^ + 0.45602 z^l^ - 2.2796 zH^ 
+ 23.231 z^l^ + 0.024306 z^ + 1.1007 zHj - 8.8726 z^l] 

[4.4496 z - 20.979 z^ + 2.5931 zl,^ - 3.0558 z^;^ 
+ 0.70052 zll - 0.49414 z^] , (44) 



^771 -1 

r<MOM ^ 

+ a, [-1.0 z - 0.5 z2 + 0.16667 z^- 0.083333 z* + 0.05 2^5 

-l.OzZ,] 
+ a^ [-7.6458 z + 0.86458 z^ - 3.5129 z^ + 13.913 z^ - 36.828 z^ 

- 6.3264 zl, + 0.10417 z^l, + 2.3866 z^Z^ + 7.3259 zH, 

- 71.347 z% - 2.0 zll + 1.0208 z'^ll + 1.1956 2^/2 

- 6.706 z*Zf_^ 30.991 z^/^j 

+ al [-59.008 z + 48.743 z^- 41.464 z/^- 5.4235 z% 

- 18.829 zll + 9.1024 z'^ll - 4.0556 zl^ + 3.5697 z^l^] , (45) 

where Iz = log{—m?/^P') and we have evaluated the coefficients in the series in 
z with Uf — 4. 

To illustrate the quality of these expansions, we have plotted (see figures [^ 
and||) the ratio of the 1, 2 and 3 loop coefficients of C^*-"^ and C^ as functions 
of 1/z = —fj^/m^ in the Landau gauge and for simplicity with Uf = 4 for all 
values of z. The circles in the plots for 1 and 2 loops correspond to the exact 
results from |4^. The convergence of the small mass (corresponding to large 
negative values of 1/z) expansions is good for 1/z < —4, where the expansions 
for higher orders of z are almost indistinguishable as well among each other as 
well as from the numbers received from the exact 2 loop propagator. On the 
other hand, due to the z ln(z)*, i = 1, . . . ,1 ( where I is the number of loops) 
terms, the MOM coefhcients are approaching the corresponding values in the 
Rl-scheme for increasing l/|z| only very slowly. This makes the Rl-schcme as 
an approximation to the MOM-scheme for the c quark useless. 
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Figure 1: The ratio of the RI and MOM scheme conversion functions C^f'^/C,^ 
as functions of —^"^/m?. Shown are the coefficients in the expansion in Os/tt as 
expansions to order z to z^ (z^ for 3 loops). Note that in Landau gauge C2 = 
, for 2 loops some numeric values for the exact mass dependence are shown as 
well. 
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Figure 2: The ratio of the RI and MOM scheme conversion functions 
^m'^^/^nii ^^ functions of —iP/w?. Shown are the coefficients in the ex- 
pansion in tts/TT as expansions to ordarcz to z^ {z^ for 3 loops). For 1 and 2 
loops some numeric values for the exact mass dependence are shown as well. 



4.2 Four Loop Quark Anomalous Dimensions 

We start from the MS scheme. The quark mass anomalous dimension was 
computed at four loops quite recently in |p3, pa| and reads 



i>0 



.(0) 



(46) 



U 16 



202 
~3~ 



Uf 



20 



(47) 



„(2) - 



64 



1249 + Uf 



2216 160 



27 



140 



(48) 



v(3) 



256 



4603055 135680 



162 27 

91723 34192 



C3-88OOC5 



Hf 



nj 



27 



18400 ^ ■ 

C3 + 880C4 + ^— Cs 



5242 800 160 ■ 

^43" + ~9"^'"~^' 



9 



332 64 

243 ^27' 



49) 



The result for the quark field anomalous dimension was found by one of 
the authors in the course of computing 7™ |5J] and read for the QCD case in 
Landau gaugeg 

72(a.)^-^7^')a:+\ 

i>0 



72 



(2) 



64 



72 



^M 


67 
3 


+ «/ 


4 
~3 




> 




20729 
36 


79 ■ 

-yC3 


+ '^Z 


55 
9 


0' 


2 


'20' 
.27. 



„(o) 



(50) 



(51) 



*72 is gauge dependent and in the Landau gauge 72"' = 0; the results for S\J{N) group in 
general covariant gauge are given in Appendix C. 



16 



(3) 
72 



256' 



2109389 565939 , 2607 ^ 

C3 + —^ C4 



162 324 

162103 2291 



761525 
1296 



81 



79 160 

^C3-yC4-— C5 



3853 160 ■ 



140 
243 



(52) 



For completeness we also give the field anomalous dimension for the case of 
QED with Uf different fermion species^: 



QED{0) 
72 



(53) 



72 



QEDil) ^ J_ 
" 16' 



[-2] 



(54) 



QED (2) ^ J_ 
'2 64 



= ^ + ^f[^\+^) 



20 



(55) 



72 



QED (3) _ 1 



256 



1027 



400(3 + 640 C5 



f-... 



+ 



n^ 



304 
~9~ 



32 C3 



280 



(56) 



In order to compute the corresponding anomalous dimensions for the RI and 
RI' schemes, one just needs to make use of Eqs. ( p6| ) to ( p9| ) in combination with 
the three loop conversion functions ( [3^j37| ). As a result for the QCD case we 
get for the RI scheme: 



7 



RI(0) _ 



(57) 



7 



RI(1) _ 



16 



126 + 71 



52 



(58) 



^For the QED case, all gauge dependence is in 72 i where ^l is defined in Appendix IaI 
(§£ = for Landau gauge) 
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% 



RI(2) ^ J_, 

64 



20911 3344 



92^ 



81 



^f 



18386 128 , 
27 9 ^ 



(59) 



vRI(3) 



256' 



300665987 15000871 , 6160 ^ ' 

Ca + ^^ C5 



648 108 
7535473 627127 



+ ni 



108 54 

670948 6416 ^ " 
■Ca 



243 



27 



4160 , ■ 

Ca + ^^ Cs 



18832 



729 



(60) 



The corresponding equations for the RI' scheme are: 



1: 



RI'(O) 



(61) 



1'. 



Rl'(l) _ 



16 



126 



Uf 



52 



(62) 



7, 



RI'(2) ^ J_. 

64 



20174 3344 



17588 128 , 

\ Ca 

27 9 ^ 



n} 



856 



(63) 



vRl'(3) 



256' 



141825253 



^f 



324 

3519059 

54 

611152 
243 



7230017 , 6160 , ' 

Ca + ^^ C5 



54 
298241 



27 



5984 

27 



4160 ^ ■ 

Ca + ^^ Cs 



16024 



729 



(64) 
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The quark field anomalous dimensions for the RI and RI' schemes are given 
in Appendix C. 

Common to all our results for the quark mass and field anomalous dimensions 
in the RI and RI' schemes (see Eqs. (p3) and Appendix C) is the absence of <^4 
even at the four-loop level. This is in striking contrast to the behavior of their 
MS counterparts and calls for an explanation. In fact, a similar absence of the 
constant C4 at four loops in the gauge beta functions recently has been neatly 
explained in [p7| . 

Unfortunately, we have not been able to extend the argument of |5j] to the 
present case. Nevertheless, we tend to agree with David Broadhurst and the 
referee in suggesting that the appearance of ^4 in Eqs. ( |49| , |5^ ) is apparently an 
artifact of the MS-scheme. 

To support this idea we demonstrate below that, in a sense, the RI/RF 
schemes are "more physical" than the MS one. To illustrate this statement, let 
us define the q dependent generalization — C^^ {as, fJ.^ / — q^), with ? = 2 or m, 
— of the conversion functions C^^ by the same Eqs. (|9|-^) but without the 
condition q^ = —fJ,^. Let us in addition define 

rfla^MV-?') = ^^'^log\c^'{as,^^y~q')\ (65) 



2 9 , 



C^'{as,f,'/-q') . (66) 



Using these definitions, the following relations hold: 

• Boundary condition: 

C'f^(a..,l) = Cf^(a,). 

• Evolution equations (they follow directly from the evolution equation of 
the fermion propagator): 

..2 Q /^Ri „, riRi .,2 ^ T^Ri n iR'7\ 

or, equivalently, 

• Relations between the Tf^ and the 'j^^^ (can be received by combining 
Eqs. (|3j6|)and(§|,||):' 

7f^(a,) = rf(a„l). (69) 

In fact, the functions T^^ happen to be both scale and scheme independent. 
In Eqs. dl9-pOh it is understood that the fermion propagator is defined in the 
MS-scheme. However, one can easily check that even if the fermion propagator 
would be defined in any other (mass-independent) scheme, then the resulting 
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change of the C^^ would amount to a rescaling by a q-independent factor which 
can not, obviously, change the very functions T^^ . 

We conclude that the T^^ are physical scheme invariant quantities, at least 
within the framework of perturbation theory. Thus, due to the relation (p9), 
the absence of (^4 in the quark mass and field anomalous dimensions in the RI 
scheme should be considered as a phenomenon which is no more puzzling than 
the similar absenceH of ^4 in the four-loop total cross-section of e+e^ annihilation 
into hadrons calculated in massless QCD |58, p9| . 

The same reasoning is fully applicable to the case of the RI' scheme. 



4.3 NNNLO relation for the RI quark mass and the RGI 



mass TTin 



It is customary to solve the RG equation (jg) for the quark running mass m.(/i) 
as follows 



to(^o) c(as(Mo))' 
where (x stands for either as(^)/7r or as(/xo)/7r) 

llm{x')' 



(70) 



c{x) 



exp 



dx' 



(3{x') 
(xr{l-f (7-i-/3i7o)a; 



(71 - /3i7o) + 72 + /3i 70 - /3i7i - /3270 

1 - 1 - - 2 - - 

g (71 - Pilof + 2 (71 ~ /3i7o)(72 + A 70 - Piii - /327o) 



+ 0{x^)) 



+ :7 ( 73 - /3i 70 + 2/3i/327o - /537o + /?! 7i - f^il\ - P1I2 



(71) 



Here 7^ = 7™ //3o, Pi = Pi/Po, (i=l72,3) and pi are the coefficients of the 
QCD beta-function as defined in Eq. (w. The four- loop beta- function recently 



has been computed in |56 with the result 



1 / 2 

Po= 4 11-3"/ 



/3^ ^ 1 ( 102 - -n, 
^ 16 V 3 ^ 



P2 
P3 



1 /2857 5033 



64 V 2 18 

1 / 149753 
256 V 6 



325 

? 

54 



3564C(3) - 



1078361 6508 



50065 6472 



162 



81 



C(3) 



162 
1093 . 



27 



C(3) 



(72) 



V 



729 



'/ 



^We do not count the well-understood 7r^-term arising due to the analytical continuation 
to the physical region of energies. 
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Eq. ( [70[ ) directly leads to the RG invariant j-L independent mass rhq 

„ _ mq{^l) 



(73) 



' c(a,(^))' 

An important property of rhq is its /i and scheme independence. The latter 
follows from the fact that rhq could be alternatively defined as follows 



lim 



M) 



(74) 



and from the well-known universality of the one loop coefficients of the quark 
mass anomalous dimension and the /3-function. 

Evaluating the four loop approximation of the c-function in the RI and 
RI' schemes, we can state our results for the conversion functions as a relation 
between the RG invariant mass rh and the masses mF^^ and m^^ . For n/ = 3, 4, 5 
Eq. ( [73| ) assumes the form: 



m 



(3) 



,RI 



1 + 7^ 
47r 



Utt 



7221 /<^ 
" 81 J Utt 

88484924345 
12754584 

3089567 



972 



Ca 



2521517 
13122 



18640 
81 



536/ 



(75) 



m(4) 



m 



RI 



(v) 



1+? 

47r 



V47ry 



17606 

' 1875 ^ 

3819632767 952 
21093750 ^15"' 

8512503162869851 
1423828125000 

1035345331 



337500 



■C3-304C5 



(76) 



m 



(5) 



tRI 



(v) 



47r 



Utt 
Utt 



15926 



1587^ 

2559841211 4696 

15111414 ^ 69 

4334826270205387 

863345304648 

3889063057 



+ 



1314036 



■C3 



26960 
69 



(77) 
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5 Conclusions 

In this paper we have analyticaUy computed the first few terms of the high- 
energy expansion of the three-loop quark propagator. These resuhs have been 
used to find the NNNLO conversion factors transforming the MS quark mass 
and the renormalized quark field to those defined in the RI scheme which is more 
suitable for lattice QCD calculations. The newly computed NNNLO corrections 
are numerically significant and should be taken into account when transforming 
the RI quark masses to the MS ones. 

We also have presented the four loop results for the quark mass and wave 
function anomalous dimensions in the RI and RI' schemes. Unlike the case of 
MS-scheme, the results display a striking absence of the C4 irrational constant 
even at four loops. This could be attributed to the fact that the RI /RI' quark 
mass and field anomalous dimensions could be defined in a scheme-invariant 



way (see subsection (L2) for more details). 

In principle, the knowledge of N^LO conversion factors would be useful to 
even better control the convergence of the perturbation series. Unfortunately, 
such a calculation requires the knowledge of the quark propagator at four loops 
— a problem certainly out of the range of present calculational techniques. 
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7 Note 

In [Q the NNNLO MS — RI conversion relations have been used to transform 
the lattice results for the RI light quark masses into those for the MS ones. In 
[ pll the results for the (QED) fermion mass and field anomalous dimensions ( 
Eqs. p4,B3, p&j and Eqs. (7-9) of [p4|) have been reproduced within an entirely 
different approach. 
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A Quark Propagators 

Below we list the full three loop results for the quark propagator^ computed in 
general covariant gauge with the tree gluon propagator 



1 



{gp.:^ - (1 - iL)qi,qu/q^)- 



For SU(A^) gauge group colour factors have the values Ca = N, Cf = 
[N"^ - 1)/(27V) and T = 1/2. The QED case is obtained with substitutions 
Ca ^ 0, Cf ^ 1 and T ^ 1. 
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B Conversion Functions 

The full gauge dependent conversion factors read 
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'Note that these results are expanded in a^ = ^^ and not in j^ as in Eqs. (|34J) to (BTJ). 
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C Anomalous Dimensions 

The quark mass and field anomalous dimensions for general SU(iV) and for the 
MS, RI and RI' schemes are listed below. See Eq. (HI) for the conventions. For 
the MS case also the field anomalous dimension 72 is given for general gauge 
and SU(iV), while for the RI and RI' only the Landau-gauge formulae are given 
(for which all 72 are 0): 
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